Phase-ordering of conserved vectorial systems with field-dependent mobility 



Federico Corberi* 

Dipartimento di Scienze Fisiche, Universita di Napoli and Istituto Nazionale di Fisica della Materia, Unitd di Napoli, Mostra 

d'Oltremare, Pad. 19, 80125 Napoli, Italy 

Claudio Castellanot 
The Abdus Salam ICTP, Strada Costiera 11, P.O. Box 586, 34100 Trieste, Italy 



The dynamics of phase-separation in conserved systems with an O(N) continuous symmetry is 
investigated in the presence of an order parameter dependent mobility M{<p} — 1 — acf> 2 . The model 
is studied analytically in the framework of the large- JV approximation and by numerical simulations 
of the JV = 2, N — 3 and JV = 4 cases in d = 2, for both critical and off-critical quenches. We 
show the existence of a new universality class for a — 1 characterized by a growth law of the typical 
length L(t) ~ t 1 '* with dynamical exponent z = 6 as opposed to the usual value z = 4 which is 
recovered for a < 1. 
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I. INTRODUCTION 



The phase-separation kinetics of conserved systems quenched from a high temperature disordered state into the 
ordered region of the phase diagram is usually modeled by the Cahn-Hilliard JjJ equation for the order parameter 
field <fi(x, t) = {4> a (x, t)}, (with a = 1, . . . , N) 



d<j> a (x,t) 
dt 



= V <^ M(»V 



6F(cf>) 



(1) 



In Eq. (0) F{4>) represents an O(N) symmetric free energy functional with ground state in (f) 2 /N = (1/N) J2a=i = 1 
and M{4>) is the mobility. Although in the original work of Cahn and Hilliard a constant M was proposed, it has 
been subsequently argued H that, for scalar systems, a field-dependent mobility 



M{(j>) oc 1 - a(T)cj) 2 



(2) 



is more appropriate; a(T) — > 1 for temperature T ~ * 0, while a(T) — * for T — > T c . For vectorial systems it is natural 
to generalize Eq. (0), considering a mobility 



M(0)txl-a(T)| 



(3) 



Phase-ordering in presence of non-constant mobility has recently been studied with JV = 1 (g|JJ], showing a richer 
behavior with respect to the usual case with constant M. From a microscopic point of view it is easy to understand the 
existence of two different coarsening mechanisms in the separation of binary mixtures, which is the typical realization of 
a scalar conserved phase-ordering system. The first is surface diffusion, namely the diffusion of molecules along domain 
walls in order to minimize the interfacial energy, and is associated with a growth law of the domain size L(t) ~ t 1 /" , 
with z = 4. The second consists in the evaporation of molecules from large curvature regions of domain interfaces 
and their subsequent diffusion through the other phase towards less curved interface portions. It is often called also 
Lifshitz-Slyozov or evaporation-condensation mechanism; the growth exponent is z = 3 @]. Since bulk diffusion is an 
activated process and surface diffusion is not, an interesting phenomenology occurs when the temperature is changed. 
For high T (but still T < T c ) one observes only bulk diffusion (z = 3), which is faster than the other mechanism. 
When the temperature is lowered, bulk diffusion is strongly suppressed due to its activated nature; it is therefore 
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possible to observe a preasymptotic regime dominated by surface diffusion [z = 4). This crossover is reproduced by 
the continuum equations with non constant mobility , when the parameter a is varied. At the coarse-grained level 
of the continuum equations, this again can be understood by observing that while for shallow quenches (a <C 1) M((f>) 
remains finite in the whole system, the situation is different when deep quenches are considered. When a — I the 
mobility vanishes within domains and suppresses bulk transport; since M(<f>) ~ I on domain boundaries, the motion 
of interfaces is instead unaffected. 

Despite the formal analogy of Eq. (|f|) for scalar and vectorial fields the underlying coarsening mechanisms are in 
principle different, as a consequence of the different symmetry of the ground state. However for N < d one has stable 
localized topological defects, of which domain walls are the scalar counterpart. With constant mobility defects play 
in vectorial systems exactly the same role of surfaces for N = 1: Phase-ordering proceeds for long times by reduction 
in their typical radius of curvature (if they are extended) or via mutual annihilation of defect-antidefect pairs (for 
point defects). This analogy is however not complete: It is not clear what are in the vectorial case the coarsening 
mechanisms corresponding to surface and bulk diffusion; it is not even clear that two different mechanisms must exist 
at all. Furthermore, when N > d no stable localized defects are present and the coarsening process has evidently a 
different nature, even if the growth exponent is the same. 

In this paper we consider the phase-ordering kinetics of a vectorial field with a non constant mobility. We report 
results obtained numerically in d = 2 for the case iV = 2, as a paradigm of systems with topological defects, and 
the cases N — 3 and 4, where stable defects are absent. We also study analytically the soluble large-iV model. 
Interestingly, we find that the global picture is the same from N = 2 up to N = oo, and perfectly analogous to the 
scalar case. More precisely, while for a -C 1 the power growth law L(t) ~ t 1 !" is obeyed asymptotically with z = 4, 
with a — 1 a new value of the dynamical exponent z = 6 is observed asymptotically, corresponding to a lower growth 
rate, in analogy with the case N = 1. For a < I the z = 6 exponent is expected to be observed preasymptotically 
before a crossover leads to z = 4. These results are found both for critical and off-critical quenches. 

From these results it turns out that the similarity between ordering in scalar and vectorial systems is very strong. 
Also for N > 1 the dynamic exponent z is changed by the depth of the quench in the presence of a non constant 
mobility. This close analogy is not limited to fields supporting topological defects but is instead valid for any number 
of components of the order parameter. 

The article is divided as follows: In Sec. 2 the model is described; in Sec. 3 the solution of the large-iV model is 
presented; in Sec. 4 the results of numerical simulations of two-dimensional systems with iV = 2, JV = 3 and N = 4 
are presented and compared to the cases N = I and N = oo; in Sec. 5 we summarize the results and draw our 
conclusions. 



II. THE MODEL 



We consider a system with a vectorial order parameter <p(x,t) initially prepared in a configuration sampled from 
an high temperature uncorrelated state with expectation values < <fr a (x,0) >= N^m a and < <j) a (x 7 0)(f>p(x' 7 0) >= 
A5 a p5(x — x'). For t > the time evolution is governed by the noiseless Langevin equation ([!]), where M(0) is given 
by Eq. (0), and F{<p} is assumed in the Ginzburg-Landau form: 



F{4>} = J dx 



ii V( /,| 2 - V + — w> 2 ) 2 

2 1 1 2 4iV v ; 



(4) 



With these positions the equation of motion for the order parameter field reads 



d(j) a (x,t) 
dt 







V 




[ N _ 





-V 2 4> a {x : i) ~ 4> a {x,t) + ^cp 2 (x,t)(j) a (x,t) 



(•5) 



The central quantity in the study of the dynamical process described above is the structure factor, namely the Fourier 
transform of the pair connected equal time correlation function, defined by: 



C a (k,t) =< cj) a {k,t)(f> a {-k,t) > -Nm 2 a S(k) 



(6) 



In (||) and hereafter < ... > means ensemble averages, that is over the stochastic initial conditions, since Eq. (||) is 
deterministic. If the m a are equal (for critical quenches) all the components of the structure factor are equivalent 
and the index a will be dropped. 
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III. THE LARGE-TV MODEL 



In this Section we present the analytical solution of the vectorial model in the limit of an infinite number of 
components of the order parameter field. Although it has been shown || that the nature of the dynamical process 
is different when N is strictly infinite, the signature of this being the multiscaling symmetry obeyed by the structure 
factor, the global picture provided by this model is often qualitatively adequate and predicts the correct dynamical 
exponent z for N > 1. Here it will be shown that the non trivial behavior of systems with non constant mobility, 
which is due to the interplay between the different coarsening mechanisms, is fully reproduced by the large- AT model. 

Let us consider a process where symmetry breaking can be induced along one direction by virtue of an asymmetric 
initial condition 

mi ^0 mp = {/3 = 2,...,N} (7) 

In the large- N limit the evolution equations for the longitudinal and transverse components of the structure factor 
(see Appendix) read 

dCi(k,t) = 2 {i- a [S ± (t)+m 2 }}k 2 [-k 2 + l-S ± (t)-3m 2 )]C 1 {k,t) (8) 



dt 

ac±{k,t) 
at 



= 2{1 - a[S±(t) + m\]}k 2 [-k 2 + 1 - S±(t) - m\)] C^{k 7 t) (9) 



where C\ refers to the correlations along the symmetry breaking direction, a. = 1, and Cj_ is the structure factor 
along one of the equivalent transverse directions. S±(t) can be computed self-consistently through 

/dV 
j^jCxiKt) (10) 
,k|<g V Z1T ) 

where q is a phenomenological ultraviolet momentum cutoff. 

Eq. (|]) together with the self-consistency relation ([l0|) governs the dynamics in the large- N model. These equations 
apply to off-critical and critical quenches (mi = 0) as well; in the latter case only one equation is required since Eqs. (||) 
and @ coincide. Notice that Eq. @ for Cj_ does not contain C\ and therefore, with an high temperature disordered 
initial condition Cj_(k, 0) = A can be formally integrated yielding 

CMM) = Ae^ 4 ^ 2 ^ (11) 



where the lengths 



\\l', = \2j {l-a[S^{T)+m 2 ]}dr^ (12) 



C(t) = hj{l- Sx(t) - m 2 }{l - a[S±(r) + m 2 }}d T ^ 



1/2 



(13) 



have been introduced. For long times the structure factor is sharply peaked around 



km = ffj (14) 



allowing a saddle point evaluation of the integral over momenta in Eq. ( |l0| ) which yields 



S(t) - A£- d (t)e^w (15) 

For long times convergence towards thermodynamic equilibrium requires the order parameter to approach the mini- 
mum of the local part of the free energy (^) which is for <j> 2 = N. Hence, letting S±(t) ~ 1 — m\ asymptotically in 
Eq. (Ha) one has 
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A(t) ~ ^ ; (16) 

[4d In 



Two different cases must then be considered, namely a < 1 or a = 1. For a < 1 one immediately obtains from Eq. ( jl^ 
A(t) ~ t 1 / 4 which yields C(t) ~ (tint) 1 / 4 due to Eq. (|13). Hence the physical length L(i) ~ ^m 1 (t) associated to 
the peak of the structure factor grows as Lit) ~ [i/ lnjft)]*/ 4 . With a < 1, therefore, the asymptotic behavior is the 
same as for a large- N system with constant mobility |Q, as expected. Notice the logarithmic correction with respect 
to the power law growth obeyed by systems with finite N which is due to the multiscaling symmetry of C(k,t); the 
dynamical exponent z = 4 is known to be correct for all physical vectorial systems. 

Let us consider now the case a = 1. In this case, by matching Eqs. (^2|), (|l^) and (|l^) one finds a different solution 
characterized by A(t) - t 1 / 6 and C(t) ~ ^(lnt) 1 / 4 , which yields L{t) ~ t 1 / 6 /^*) 1 / 4 . With o = 1, therefore the 
vanishing of the mobility in equilibrated regions slows down the dynamics and changes z from 4 to 6, similarly to 
the scalar case where one goes from the Lifshitz-Slyozov evaporation condensation mechanism, associated to z = 3, 
to z = 4. It must be noticed however that, apart from this analogy, the physics of the coarsening process is very 
different due to the absence of stable topological defects for N > d and there is no clear indication on the nature of 
the transport mechanism associated with z = 6. 



IV. SYSTEMS WITH FINITE TV 

In this section we present the results of the numerical solution of Eq. (|^). Our aim is to present a rather complete 
description of the effect of a non-constant mobility in systems with different number of components and to compare 
the results with what is known for N = 1 and with the above discussed large- TV model. In order to fulfill this program 
we have chosen d — 2 and TV — 2 or TV = 3 and 4 so that we can analyze systems with and without stable topological 
defects. 

The numerical solution of Eq. (|^) is obtained by simple iteration of the discretized equation on a 512x512 mesh, 
except for the off-critical N = 2, a = 1 case which is computed on a 256x256 lattice. Each quench is averaged over 
three different realizations of the initial conditions. The numerical solution of Eq. (||) is particularly delicate for a < 1, 
because in this case a small error in the calculation of cf> in the bulk may give rise to a negative M(<f>), thus producing 
a spurious instability. We have avoided this problem by letting M{<f>) = |1 — a(f> 2 \, and tested that this does not 
produce any significative difference in the results. 

The characteristic growing length L(t) is obtained from the structure factor as Lit) — fc^ 1 (t), where 

JdkkC(k,t) 

* l( ' ) = JdkC{k,t) (17) 

is the first moment of C{k,t). 



A. TV = 3 and TV = 4 

In Fig. (0) the behavior of L(t) is compared in the two cases a = and a = 1, for critical quenches with N = 3 
and N = 4. After the linear instability is over (which corresponds to times up to ln(t) ~ 4), the system enters the 
asymptotic regime which is characterized by a power law growth of L{t). For a = one has z — 4 with good accuracy 
almost immediately, while for a — 1 a convergence towards z — 6 is observed. Best fit estimates for times larger than 
t = 3000 yield z = 3.92 and z — 3.89 for a — and N = 3 and N — 4 respectively; for a = 1, instead, one has 
z = 5.68 and z = 5.81 in the corresponding cases. We notice, by the way, the remarkable superposition of the curves 
with N = 3 and N = 4. 

It is interesting to observe the evolution of the structure factor, which is shown in Figs. (Q) and (BH, for N = 3, 
at different times for a = and a = 1, respectively. Here one observes, immediately after the linear instability 
corresponding to the exponential growth of a peak at constant wave vector, the formation of a rather well developed 
power law decay of C(fc,t) for large fc; this tail is initially formed at very large values of k 7 subsequently it moves 
towards lower wave vectors and then gradually disappears as time increases, being replaced by a faster (exponential) 
decay of the structure factor. As Fig. (||) shows the exponent of this decay is very close to 5 in accordance with the 
existence of a generalized Porod's tail, i.e. C(k,t) ~ ^-(N+d) ^ j g we vj known that the Porod's tail is associated 
with the presence of localized topological defects. The presence of such a tail in Figs. (Q) and (||), suggests therefore 
that defects are formed at early time and eventually decay so that the power law disappears asymptotically. This 
phenomenon is displayed, with similar features, both for a = and for a = 1, and will be studied elsewhere M. 
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In Figs. (0) and (||) the data collapse for C(fc,t) are presented. Here 5 curves corresponding to about 1.5 decades 
in time are superimposed by plotting L(t)~ d C(k, t) against x = kL(t). We observe that for a — the data collapse is 
very good for x < 2, whereas it becomes progressively less accurate for increasing values of x. This effect is probably 
due to the remnant of the power law tail that, for the times covered by the numerical solution, has not completely 
disappeared, as can be seen from Fig. (||). In any case the data are consistent with the existence of the scaling 
property both for a = and a = 1. In Fig. M) the scaling functions of the cases a = and a = 1 are compared. The 
scaling function is obtained here as L(tu) C(k,tM), tu being the longest time of our computation. It has been 
shown in Q that no significative differences are observed in the scalar case between the two scaling functions, a fact 
that indicates the independence of the morphology of the growing phases with respect to a change of the coarsening 
mechanism from bulk to surface diffusion. Fig. (J6|) shows that the same property is shared by our vectorial model. 
Notice also that C(k,t) ~ k 4 for small k pTo[ ] . 

B. TV = 2 

We present in the following the results of the numerical simulation of the case N = 2. In Fig. (|t]) the behavior of 
L(t) is plotted versus t in the two cases a = and a = 1. From this figure we conclude that a power law growth with 
z = 4 or z = 6 for a = and a — 1 respectively, fits very well the data. Linear regression analysis for t > 3000 yields 
z = 4.20 and z = 6.26 for a = and a = 1, respectively. 

The structure factors at different times for a — and a = 1 are shown in Figs. (||) and ([}]), respectively. Here one 
observes again the formation of the Porod's tail but, differently from the cases N = 3 and N = 4, this pattern is 
maintained asymptotically, a fact that reflects the stability of the topological defects in the late stage of the dynamics. 
In Fig. ( [To|) and ( pd|) the scaling plots for C(k,t) are presented. No compelling evidence can be obtained by the 
analysis of this figures about the possible existence of a scaling breakdown, which has been suggested with constant 
mobility for d = 2 and N — 2 nevertheless one observes a worse data collapse with respect to the case N = 3 
and 4 for large x, both for a = and a = 1. For x < 2, on the other hand, the collapse is rather good. In Fig. ( [l2] ) a 
comparison is presented between the quantity L(tM)~ d C(k,tM) (which, if scaling exists, corresponds to the scaling 
function) in the cases a = and a = 1. Here one observes a superposition of the two curves for x < 2. For x > 2, 
however, they become more and more distant with x. The situation is different from the cases N = 1 and N > 2 
where the scaling functions practically coincide in the whole range of x values. From the observation of Figs. (|iol) , (|Tl|) 
and (jl|) we argue that, for x < 2 the N — 2 model behaves very similarly to the cases N = 1 and N > 2, since we 
find a superposition of curves at different times in the scaling plots (|ici|), jll|), and very similar "scaling function" for 
a = and a = 1. For x > 2, however, one observes a much worse superposition in both cases. It is not presently clear 
if the origin of this difference is due to preasymptotic effects in the simulations or to scaling violations. 

In Fig. ( |l3| ) the effect of a finite m a is considered and the growth of the characteristic length L(t) is shown vs. t for 
a = 1 both for critical and off-critical quenches. The asymmetric quench is realized by taking m a = 0.65/V2, Va. The 
solution for the off-critical case has been obtained on a 256x256 lattice in order to speed up the computation since, 
as the figure shows, much longer times than in the critical case must be reached in order to observe the asymptotic 
behavior. For sufficiently long times one observes z = 6 in both cases, suggesting that the same coarsening mechanism 
is at work both for critical and off-critical quenches. 

V. CONCLUSIONS 

In this paper we have considered the effect of an order parameter dependent mobility on the phase-ordering process 
of a vectorial system quenched in the ordered region of the phase diagram. This problem has been studied both 
numerically, by considering the cases N = 2, N = 3 and N = 4, and analytically in the large- TV limit, for critical and 
off-critical quenches. We have shown that for a = 1 the dynamics is slowed down with respect to the case a = 0, due 
to the vanishing of the mobility in pure phases. A power law growth of the typical length L(t) ~ t x / z is still obeyed 
in this case but with z — 6. This behavior is exhibited in vectorial systems for all the values of N considered and 
does not depend on the symmetry of the quench, nor on the presence of topological defects. We conclude that this 
case falls within a different universality class with respect to a = 0. For intermediate values of a the system behaves 
as if a = asymptotically but a preasymptotic growth with z = 6 is expected for a sufficiently close to 1. 

These features are very reminiscent of those of scalar systems, where the dynamical exponent changes from z = 3 
to z = 4 when a = 1. It is however important to point out that the scalar models studied in Ref. || and [Q are not 
perfectly equivalent. The difference is in the form of the free energy functional which is of the Ginzburg-Landau type 
(quartic polynomial) in Ref. 0] and of logarithmic type in B . While the polynomial can be seen as the expansion of 
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the logarithmic form for T not too far from T c , for lower temperatures only the logarithmic form provides a correct 
coarse-grained description of the microscopic dynamics. This is confirmed by what happens for off-critical quenches 
to T = 0. When the quench is sufficiently asymmetric so that the minority phase forms a non percolating pattern 
of isolated droplets, surface diffusion alone cannot drive the phase-separation to completion. The system therefore 
remains pinned in a configuration out of equilibrium. This pinning is correctly reproduced by the continuum equation 
with the logarithmic free energy fll2| ; no such pinning is found with the quartic polynomial |fT3| . In this last case the 
deep off-critical condition has only the effect of changing bulk-diffusion to a slower mechanism (z — 4) which has been 
referred to as bulk sub-diffusion. The dynamical exponent for bulk sub-diffusion turns out to be the same as that for 
surface diffusion: We do not know whether this equality is accidental or has a deeper meaning. 

In our study we have considered the generalization to a vectorial order parameter of the model with Ginzburg- 
Landau free energy. Also in this case we find two universality classes for a — 1 and a < 1. This is perfectly analogous 
to what happens in the scalar case. Moreover the scaling functions for a = and a — 1 are identical. This leads us 
to conjecture that a mechanism of the same nature drives phase-ordering for all values of a. The only special feature 
of the a — 1 case is that the mechanism is slowed down, exactly as for N — 1 one goes from bulk diffusion to bulk 
sub-diffusion, but the nature of the mechanism is still of bulk type. 

It would be interesting to check whether a sort of generalized surface diffusion exists for N > 1. This could be 
detected by solving the equation of motion for a = 1 with a logarithmic free energy for a vectorial order parameter. An 
indication that this could be the case is provided by the analytical solution of the large- N model with a logarithmic 
form of the free energy [Q . 
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APPENDIX 

Defining the fluctuation field ip(x,t) as 

(f>i(x,t) = N^mx +ip(x,t) (18) 
and inserting into Eq. (^|) one obtains the pair of equations 

d[Nh mi + i,(x, t)] = v I {1 _ o[ 1 £ ^ (a . t)+mj + ^ t) + 2N ^ mi ^ t)]} 



dt I N 

(3=2 



V{-V 2 V(z, t) -N^mi- ip(x, t) + —[N^rm + ^{x, t)} ^ <f%(x, t) 



13=2 

i 3 1 1 

N?m\ + im\ip{x, t) H -m^ix, t) + —ip 3 (x, t)} \ (19) 

f 1 N 

V I {1 - a[— <t>U x > *) + m i + ^O". *) + 2N^-mnP{x, t)}}W{-W 2 c/)(}(x, t) 

[ 0=2 

1 N 2 1 ] 

-cj>p{x,t) + —J2<Pl{ x ^)M x ^) + \ m \ + —rmi4>{x,t) + —^ 2 {x,t)}(f,p(x,t) \ (20) 

a=2 ^ 2 ) 

with p^l. 

In the large- N limit, summing over vector components averages the system over an ensemble of configurations and 
hence 



dt 



I i N 

lim -|^(aj,t)| 2 = lim - ^ ^(x,t) =< <j%(x,t) >= S ± (t) (21) 



13=2 



G 



where translational invariance has been assumed and S±(t) does not depend on /3 due to the internal symmetry. To 
leading order in N one obtains: 

= {1 - a[S±(t) + m 2 ]}V 2 {[-V 2 -1 + S ± (t) + 3m^(a;,t)} (22) 



dijj(x, t) 
di 

dt 



= {1 - a[S±(t) + m\]}k 2 {[-V 2 - 1 + S±(t) + m\](t>p{x, t)} (23) 

Introducing the longitudinal and transverse part of the structure factor, namely C\(k, t) =< ip(k 7 t)ij){—k, t) > and 
C±(k,t) —< <t>p(k,t)<t>p(—k,t) >, which is independent on (3 due to internal symmetry and Fourier transforming one 
..l.i ;., , i... !■, .11 — iations 

= 2{l-a[S^(t)+m 2 ]}k 2 [-k 2 + l-S^(t)-3m 2 ]C 1 {k,t) (24) 



obtains the following pair of equation 
dd(k,t) 



dC ^ l) = 2 {1 - a [S±(t) + mj]}k 2 [-k 2 + 1 - S±(t) - m\]C ± {k,t) (25) 
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FIG. 1. The growth law of the characteristic length L(t) is shown for systems with N = 3 and iV = 4. Straight lines 
represent power laws t 1 ^ with z = 4 and z = 6. 

FIG. 2. The evolution of C(k,t) is shown at different times for a — and N = 3. The straight line represents a power law 
decay fc~ 5 . 

FIG. 3. The evolution of C(k,t) is shown at different times for a — 1 and iV = 3. The straight line represents a power law 
decay fc -5 . 

FIG. 4. Data collapse (scaling plot) for a = and N = 3. Here L(t)~ d C(k,t) is plotted against x = kL(t) for different 
times. 

FIG. 5. Data collapse (scaling plot) for a = 1 and N = 3. Here L{t)~ d C{k,t) is plotted against x = kL{t) for different 
times. 
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FIG. 6. Comparison between the scaling functions for N = 3 and a = or o = 1. The straight line represents the power 
law fc 4 . 

FIG. 7. The growth law of the characteristic length L(t) is plotted against time for a system with N = 2. Straight lines 
represent power laws t 1 ^ with z = 4 and 2 = 6. 

FIG. 8. The evolution of C(k,t) is shown at different times for a = and N = 2. The straight line represents a power law 
decay A; -4 . 

FIG. 9. The evolution of C(k,t) is shown at different times for a = 1 and N — 2. The straight line represents a power law 
decay fc~ 4 . 

FIG. 10. Data collapse (scaling plot) for a = and N — 2. Here L(t)~ d C(k,t) is plotted against x = fcL(i) for different 
times. 

FIG. 11. Data collapse (scaling plot) for a = 1 and N — 2. Here L{t)~ d C{k,t) is plotted against a; = kL(t) for different 
times. 

FIG. 12. Comparison between the "scaling functions" for iV = 2 and a = or a = 1. The straight line represents the power 
law k 4 . 

FIG. 13. The growth law of the characteristic length L(t) is plotted vs. t for a system with N — 2, for a — 1, for critical 
and off critical quenches. The straight line represents the power law t 1 ^ with z = 6. 
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